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We develop a point-scattering approach to the plane-wave optical transmission of subwavelength 
metal hole arrays. We present a real space description instead of the more conventional reciprocal 
space description; this naturally produces interfering resonant features in the transmission spectra 
and makes explicit the tensorial properties of the transmission matrix. We give transmission spectra 
simulations for both square and hexagonal arrays; these can be evaluated at arbitrary angles and 
polarizations. 



Experiments have revealed the crucial role played by 
surface wave excitations (often called surface plasmons) 
in the case of extraordinary transmission features of 
nanoperforated metallic films 1 . This experimental work 
has generated an important theoretical literature which 
can be grossly separated into two parts. A first cate- 
gory is based on ab initio fully numerical simulations 
of the scattering amplitudes of these nanoarrays; in 
this case, the interpretation of simulated spectra is of- 
ten difficulliSiMi. A second category is accordingly de- 
voted to the search for physical understanding of the 
phenomenonSiS. The main criticism that can be ad- 
dressed to these interpretative papers is that they rely 
on the a priori definition of resonances as surface modes 
at a smooth interface (i.e. no holes) which are coupled 
to the incident field via momentum matching by the hole 
array. 

The aim of this paper is to fill the intermediate gap 
with a simple model which allows both a clear physical 
description and straightforward simulations of transmis- 
sion spectra. The main attractive features of our model 
are related to the fact that it is based on a real space (in- 
stead of a reciprocal space) surface wave scattering anal- 
ysis. We stress that our work naturally produces trans- 
mission resonances and that it does so without invoking 
a priori wavevector matching between surface modes on 
a smooth interface and grating momenta, i.e. recipro- 
cal lattice vectors. Most importantly, it reveals interfer- 
ence effects between neighbouring resonances, contrary 
to the usual reciprocal space (i.e. Fourier) description 
of wavevector matching. Our description makes explicit, 
among others, the polarization dependence and the ten- 
sorial properties of the transmission matrix. It describes 
the influence of the incident angle tuning on transmission 
spectra and, eventually, it naturally relates correspond- 
ing band structures to symmetries of the reciprocal lattice 
of the array. All these aspects are, for the sake of demon- 
stration, best clarified in a rather simplified framework, 
as we will discuss below. This implies that we do not aim 
at quantitative agreement with experimental transmis- 
sion spectra. Nevertheless, we feel that our model yields 
original physical insights into the dynamics of transmis- 
sion through nanohole arrays. 

Our approach is rooted in the role played by surface 
modes on the transmission process. The transmission 
process of a plane wave through sub-wavelength metallic 



hole arrays will be cast into the context of a Huygens- 
type principle, where the array is discretized as a lat- 
tice of holes acting as point-scatterers that scatter the 
incident radiation coherently into two-dimensional sec- 
ondary wavelets. Specifically, the array is contained in 
the (x,y) plane of a cartesian (x,y,z) frame and is illu- 
minated by a paraxial plane wave of wavevector ki n and 
polarization vector ui n . This wave is described by far- 
field angles k;„/ |ki n | ~ (8 x ,9 y ,l) and two electric field 
components E ~ (E x , Ej,), decomposed into the basis of 
the TE and TM polarizations. Retaining only the th 
diffraction order, the transmission can be formulated as 
E out [A, 6] = t • Ej n [A, 6] with a 2 x 2 transmission matrix 
t[A, 9], A being the wavelength and = (0 x ,0 y ) the far 
field angles of the incoming and outgoing field. 

As schematically shown in fig.Q, we will distinguish 
two contributions to this transmission matrix 5,9 . A first 
one, tg , , wavelength dependent and proportional to 
the identity matrix, corresponds to a transmission of the 
incoming field directly through the holes, i.e. to a Bcthc- 
type diffraction regime^. A second contribution, t tt , 
corresponds to the resonant part of the transmission ma- 
trix, on which we focus hereafter. Fig.Q describes this 
resonant transmission process as a three-step process: (i) 
the incident plane wave is converted into a surface wave 
at a given point scatterer, (ii) the surface wave propa- 
gates on the surface of the array and (iii) is eventually 
re-emitted as a plane wave through the array. Then, in 
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FIG. 1: Schematic representation of the scattering process, 
including the direct scattering contribution defined at a cho- 
sen origin of the lattice. 



the spirit of Fresnel diffraction, the global resonant ma- 
trix t , evaluated at the center (0) of the array, is 
based on the near-field of the array, as the coherent re- 
summation over all the secondary surface wavelets emit- 
ted by each of the holes, acting as independent uncoupled 
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holes with elementary scattering matrices r 
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(0) = ^z(0,r J ) 

{r 3 } 



(1) 



The summation is defined on the lattice distribution of 
holes {rj}. The singular term {rj = 0} is related to 
the direct transmission channel of t , and is therefore 
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naturally excluded from this resonant part. As a periodic 
lattice, the array is coordinated by its two primitive vec- 
tors (ai,a 2 ). We restrict this paper to |ai| = |a 2 | = ao 
square and hexagonal lattices, as these symmetries are 
most popular experimentallyiii£. The position of each 
scattering hole on each of these lattices is defined as 
Yj = nai + ma.2 with (n, m) integers. The summation on 
Eq.JTJ therefore corresponds to a double (n,m) ^ (0,0) 
summation over the lattice. 

In this work, we calculate the r matrix from the dy- 
namics of a scattered surface wave characterized by a 
complex transverse wavevector |k| = (771 + ir)2)2-K /A. On 
a smooth metal-dielectric interface (i.e. no holes), with 
dielectric functions £i,£ 2 respectively, theory predicts^ 



m + ii?2 



£l£2 

£1 + £2 



(2) 



Note that this expression holds for a metal-dielectric in- 
terface, from the visible to the microwave domain, but 
not for a dielectric-dielectric interface. In practice, the 
observed transmission resonances in metallic nanoholc 
arrays are much broader than and red-shifted by typi- 
cally a few percent from the mode dispersion on a smooth 
interface&iii^. These discrepancies are most likely re- 
lated to direct Bethe-type transmission channel and to 
radiative losses of the surface waves when they scatter 
on the holes of the actual structure. In our simulations, 
we use realistic values for iji and 772 (see below). 

In our model, the polarization of each surface wave is 
taken along its propagation direction with a unitary po- 
larization vector Uj = rj/\rj\. Surface polarization plays 
a crucial role, as it determines both the incoupling effi- 
ciency Uj • Ui n between the free-space incident photons 
and the excited surface wave, as being proportional to 
their electric field overlap, and the polarization of the 
emitted radiation. The incoupling factor Uj ■ U; n = cos (p 
corresponds to a two-dimensional dipole radiation pat- 
tern for the surface wave emitted at the hole; this pat- 
tern has been experimentally observed^. The full polar- 
ization behaviour is contained in the tensorial Uj ® Uj 
nature of the elementary point-scattering matrix (<g> de- 
notes a tensorial product). 

We assume that the elementary scattering matrix r is 
spherical in the far-field so that it reads as 



r(0,r,) = /(|k|)- 



i|k||r 3 -| 
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The scattering amplitude is /(|k|) = 
s(|k|)e _I7r / 4 A/Re(|k|)/27r in the far-field approxima- 



tion |k||rj| 3> 1 of the Huygens phase, being certainly 
satisfied for A <C 27rao?7i. For simplicity, we neglect the 
frequency dependence of the shape factor s(|k|) and, in 
the point-scattering limit, replace it by a constant s. 
This limits the discussion to spatial symmetries of the 
lattice. Point-group issues matter when specific shapes 
of the scatterers are introduced. Then, the scattering 
amplitude has to include a true s(k) shape factor. The 
coherence of the surface scattering is insured by the 
fixed phase relation e in ' rj with the incident field ki n . 
Eventually, with an input plane wave Ej n = Ej n u; n and 
no polarization analysis in transmission, we evaluate the 
resonant part of the intensity transmission coefficient 
T=|| scatt ;u in | 2 . 

In our simulations, wc normalize the wavelength by 
Ao = aotji and the transmission matrix by the constant 
scattering amplitude factor s. Internal damping and 
radiative losses of a surface wave propagating on the 
lattice are quantified with a single effective parameter 
•q = 772/771- This r\ determines the convergence speed 
of the scattering summation. For simplicity, we neglect 
the frequency dependence of both r\\ and 772 and choose 
77 ~ 0.02 which correspond to a constant mean free path 
of ao/(47T7/) ~ 3 /jm for ao = 0.7 /zmit Convergence is 
then easily reached with a lattice of 80 x 80 points. If 
the damping is lower, a larger lattice should be chosen. 
The frequency dependence of t g ^ is determined by 
running simulations for 800 wavelengths positioned on 
a regular grid ranging from X/X = 0.4 to A/A = 1.2. 
This interval is very large (factor of 3) but is neverthe- 
less consistent with our model, avoiding both shorter 
wavelength where the direct transmission channel takes 
over the resonant one, and longer wavelength for which 
the far-field approximation might break down. These 
nominal parameters given, the evaluation of an element 
of the transmission matrix, at a given incidence angle, 
takes about 5 seconds for the whole spectrum. 

As a first working example of our method, we present 
in fig.® a simulated spectrum of the transmission 
coefficient of a square array over the wavelength range 
mentioned above, under plane-wave illumination at 
normal incidence 8 = 0. In this case, both TE and 
TM polarizations are degenerate and the transmission 
matrix is simply proportional to the identity matrix. 
Our model immediately shows the resonant behaviour 
of the transmission spectrum. We stress again that we 
make here no use of the usual momentum conservation 
argument specifying resonances in terms of reciprocal 
lattice; the resonances arise naturally from the summa- 
tion of Eq.JTJ). 

The position and the relative strength of the peaks 
shown in the spectrum of fig. I© can be understood from 
a Fresnel-zone type argument^. A resonance builds 
up each time the Huygens phase satisfies an in-phase 
condition, related to the definition of the reciprocal lat- 
tice. Resonances are thus expected at approximately the 
same positions predicted by the usual ad hoc dispersion 
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relation for surface waves on periodic arrays, that is at 

Arcs _ 1 
Ao V n 2 + m 2 



(4) 



for a square array and A rcs /Ao = v / 3/2v / n 2 + nm + rn 2 
for a hexagonal array&. We adopt hereafter the conven- 
tion of indexing a surface mode as [n, m]. The strength of 
the peaks scales with the corresponding inverse distance 
1 / \fn 2 + m 2 . The increase in peak transmission for the 
[2,1] mode as compared to the [2,0] mode results from 
an increase in degeneracy from 4— fold to 8— fold, lead- 
ing to a factor 2 increase in amplitude and a factor 4 in 
intensity. The spectrum in fig.© is shown in arbitrary 
units (a.u.) where, from Eq.©, one hole at one lattice 
spacing contributes one arbitrary unit. 

We note that the two-dimensional summation on a dis- 
crete two-dimensional lattice (see Eq.©) does not pro- 
vide all in-phase conditions. For the square array, for in- 
stance, the in-phase conditions derived from Eq.© read 
as 



r Ap 

Arcs 



1~Kp 



(5) 



with p a given positive integer. However, n 2 + m 2 can 
not provide all such p integers. For this reason, there 
is no resonance between the [1,1] and [2,0] modes of a 
square array, though such an in-phase condition could in 
principle be satisfied in between. This lack of resonances 
corresponds to a modification of the simple Fresnel-zone 
argument: with a two-dimensional lattice, the only mul- 
tiplicities corresponding to resonant conditions are those 
matching the discretization of the lattice. 

This picture is useful to understand the positions and 
strengths of the resonances. But we would like to stress 
that our real space-based model goes beyond that, pro- 
ducing interferences between surface modes that have a 
strong influence on the transmission spectrum. Peaks 
present asymmetries due to the interference of overlap- 
ping tails of adjacent resonances. These tails can even 




FIG. 2: Transmission coefficient of a square array at normal 
incidence, as a function of normalized input wavelength. 

interfere quite destructively as in between modes [2, 1] 
and [2,0] (see fig.©)- The interferences between the 
modal amplitudes make the actual position of the reso- 
nance A rcs /Ao slightly different from the "pure" location 
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We have also investigated non-normal incidences that 
correspond to tilts along the y axis, for both square and 
hexagonal arrays (oriented with one lattice vector in the 
x— direction). Rotation angles ranging from 0° to ~ 11.5° 
in steps of ~ 2.3° were studied. TE and TM polariza- 
tions, now distinct, are displayed in separate figures. A 
central result of this paper, shown in figs.©, is that our 
model immediately and naturally produces spectral band 
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FIG. 3: Transmission spectra for a square array (first row) 
for a)TE and b)TM polarizations and for a hexagonal array 
(second row) for c)TE and d)TM polarizations. These spec- 
tra have been calculated as a function of normalized input 
wavelength and angle of incidence when the tilt is performed 
along the y— axis - see insets. For clarity, we show a globally 
shifted curve for each angular increment of ~ 2.3°. 

structures. 

For a square array, illuminated with a TE polarized 
wave, the [1,0] resonance remains practically stationary, 
whereas it splits into a doublet for a TM polarization as 
the angle of incidence is increased. For the [1, 1] mode, 
one observes a doublet splitting for both polarization. 
For a hexagonal array, the principal resonance splits into 
a triplet with a central peak for TE polarization, whereas 
it evolves into a doublet for TM polarization. Split res- 
onances will eventually overlap, through their progres- 
sive shifts. This can lead to a destructive interference, 
as seen at an ~ 11.5° incidence angle on fig.© a)) be- 
tween the [1, 1] and [2, 0] resonances for TE polarization 
or to a constructive interference, as seen on fig.© b)) 
between the [1,0] and [1, 1] modes for TM polarization. 
One should also note that for both square and hexago- 
nal arrays, central peaks are in fact only approximately 
stationary as the tilting angle is varied. With an increas- 
ing angle of incidence, a positive transverse component 
of the wavevector k; n emerges and leads to a slight blue 
shift of central peaks. 

It is interesting to compare these simulations, and in 
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particular the splitting of resonances, with what could 
be expected from rather intuitive arguments. Intuitively, 
a nanohole array is simply characterized by the spatial 
symmetries of its direct lattice. The direction of propaga- 
tion of a given surface mode will be naturally associated 
to a particular axis of this direct lattice, with an incou- 
pling efficiency argument as presented above and based 
on the projection factor between such a propagation axis 
and the incident polarization. Thus, one forbids modes 
which propagate along directions perpendicular to the in- 
cident polarization to be excited. Within this intuitive 
frame, one predicts splitting of resonances by merely per- 
forming symmetry operations on the direct lattice as the 
array is tilted. It turns out however that following this 
line of reasoning leads to wrong predictions. 

On the contrary, it is important to realize that the 
band structures of figs.® can only be fully inferred if one 
realizes that Eq. Q , starting indeed from point-scatterers 
on the direct lattice of the array, provides resulting waves 
that are directional and that propagate along axes of the 
reciprocal lattice rather than the direct one. Therefore, 
these band structures are only consistent with the sym- 
metry arguments performed on the reciprocal lattice^, 
together with a polarization "selection-rule" killing sur- 
face waves that are allowed by symmetry to propagate 
on reciprocal lattice axes perpendicular to the incident 
polarization. 

The disagreement with the intuitive approach men- 
tioned above is not present for a square array since in this 
case, the direct lattice and the reciprocal lattice are iden- 
tical. It is however easily visible for a hexagonal array for 
which direct and reciprocal lattices do not coincide. In 
this case, if one concentrates the symmetry arguments on 
the direct lattice, the expected splitting of resonances will 
not agree with simulations. For instance, symmetry ar- 
guments applied to a direct hexagonal lattice tilted along 
its y— axis lead to the prediction of quadruplet splitting 
with no stationary mode for the principal resonance, for 
both polarizations. The same argument applied to the 
reciprocal lattice, together with accounting for the pro- 
jection factor between lattice axes and polarization, leads 
to a triplet, respectively a doublet, for TE, respectively 
TM, polarization; this is indeed observed in fig.® c)) 
and fig.® d)). 

In conclusion, we have presented a simple and straight- 
forward model based on a Huy gens- type principle. This 
work provides physical insight into the surface-plasmon- 
assisted transmission process through metal nanohole ar- 



rays, emphasizing symmetry and tensorial properties of 
the transmission amplitude. It yields band structures 
consistent with directional collective surface excitations 
propagating on the reciprocal lattice of the nanohole ar- 
ray and does not rely on any ad hoc momentum matching 
argument. 

In order to stress the core characteristics of our Huy- 
gens description, we have used a simplified model that 
addresses the resonant contribution t to the trans- 
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mission through a single interface only and that does not 
account for the direct transmission channel contribution 
-b the' ^* ur P om t- sca ttering approach amounts to con- 
sidering the holes in the metallic film to be in the far 
subwavelength limit, i.e. basically in the limit of holes 
of zero diameter. To quantitatively compare with exper- 
imental data, the following extensions are needed: both 
the dielectric constant of the metal and the scattering 
amplitude s(k) require a frequency dependence. Further- 
more, related in particular to hole size effects, this scat- 
tering amplitude should be specified and the direct con- 
tribution t to the transmission should be included, 
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providing resonance line shapes with red-shifts and red- 
taik 9 . 

Of course, a qualitative agreement between our simu- 
lations and angle dependent transmission measurements 
should be easier to obtain, in particular as far as split- 
tings of resonances are concerned. But discrepancies are 
already found at this less ambitious level when refer- 
ring to experiments performed with standardly designed 
nanohole arrays 67 . In fact, such nanohole arrays cor- 
respond to optical systems with two-interfaces, possibly 
identical in the case of free-standing films. There, an- 
gle dependent couplings between surface modes defined 
on each interface induce perturbations too strong to al- 
low even a qualitative matching with the mere single- 
interface band structure provided by our model. Nev- 
ertheless, if one could design an effective single-interface 
nanoperforated film, such a qualitative check should be 
possible. A metallic film, with very small holes, deposited 
on a given dielectric substrate with a thick titanium (Ti) 
bonding layer is likely to be a relevant candidate^. Strong 
absorbtion in this Ti layer prevents indeed any surface 
mode from being excited on this metal-dielectric inter- 
face, thus keeping only a single air-metal interface into 
play. In this framework therefore, it will be interesting 
to confront experimental transmission spectra with the 
Huygens description formulated in our paper. 
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